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Nuclear parton distribution functions (NPDFs) are investigated by analyzing the 
data on structure functions and Drell-Yan cross sections o"^'^- An important 
point of this analysis is to show uncertainties of the NPDFs by the Hessian method. 
The analysis indicates that the uncertainties are large for antiquark distributions 
at a: > 0.2 and gluon distributions in the whole x region. We also discuss a nuclear 
effect on the NuTeV sin^ 9w anomaly as an application. 



1 Introduction 

We report an analysis on nuclear parton distribution functions (NPDFs) and its 
application to the NuTeV sin^ 9w issue. Although the PDFs in the nucleon are 
now known rather accurately, their nuclear modifications are not well determined. 
It is important to determine the NPDFs not only for establishing QCD in nuclei 
but also for applications to heavy-ion physics and neutrino reactions. 

We proposed the optimum NPDFs in 2001 by a analysis of nuclear data on 
deep inelastic lepton scattering ^1^. In recent years, it becomes important to show 
reliability of the obtained PDFs. Actually, uncertainties of the PDFs have been 
investigated in unpolarized and polarized PDFs in the nucleon. Since there was 
no serious error estimate of the NPDFs, we investigated the uncertainties by the 
Hessian method in the recent version In addition, Drell-Yan and HERMES 
data are added into the data set, and charm-quark distributions are included in the 
new analysis. 

The NuTeV collaboration reported that their neutrino scattering experiments 
indicate a significant deviation for the weak mixing angle sin^ 9^ from other mea- 
surements. Because the target is the iron nucleus, nuclear effects need to be in- 
vestigated before discussing any new physics mechanisms. In particular, we study 
effects due to the difference between and dv nuclear modifications '3'4' . 

This paper consists of the following sections. In Sec. 12 the analysis method 
is explained for obtaining the optimum NPDFs. Results on the NPDFs are shown 
in Sec. |21 and their effects are discussed on the determination of sin^ Ow in Sec. 01 



2 Analysis method 

The NPDFs are expressed in terms of parameters, which are determined by an 
analysis of the data on nuclear structure functions and Drell-Yan processes. 
The NPDFs are taken as a nucleonic PDF multiplied by a weight function Wi 
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which indicates nuclear modifications: 

f,^ix,Ql)^w,{x,A,Z)Mx,Ql), (1) 

where A and Z are mass number and atomic number, and i denotes a parton- 
distribution type. The function Wi is expressed by parameters: 

, . . / 1 \ a^{A, Z) + hx + c^x^ + djX^ 

The A dependence is motivated by a simple picture of nuclear volume and surface 
contributions to cross sections. The cubic functional form for the x dependence is 
motived by the shape of F^/F^ data, and the factor 1/(1 — x)^' is introduced to 
reproduce the Fermi- motion part. 

For the NPDFs (Z/'), we take u^, d^, q^, and g"^, and the initial point is chosen 
(5q=1 GeV^. These distributions are evolved to experimental points to calculate 

= J^jiRj"*" ~ Rf''°f/{(^'j"'^''Y, where ap^" is an experimental error and Rj 
indicates the ratios F^ jF^ and (y^jyyl^^DY- Leading-order expressions are used in 
the theoretical calculations. By the analysis, we obtain the optimum NPDFs 
and a Hessian matrix iJ. Using this matrix, we calculate the uncertainty of the 
distribution j'^ix) by 




where is a parameter, and ^ indicates the parameter set for the minimum 
The Ax^ determines a confidence region, and it is taken Ax^ = 10.427 for nine 
parameters "T. It corresponds to the one-cr-error range. 

3 Results on nuclear PDFs 

The experimental data, which are used for the analysis, are taken for the nuclei: 
deuteron (D), helium-4 ('^He), lithium (Li), beryllium (Be), carbon (C), nitrogen 
(N), aluminum (Al), calcium (Ca), iron (Fe), copper (Cu), krypton (Kr), silver 
(Ag), tin (Sn), xenon (Xe), tungsten (W), gold (Au), and lead (Pb). The total 
number of the data is 951. The minimum ^ becomes Xmin^l^^Q-^ i'^ the current 
analysis. 

Among many nuclear data sets, we show only two examples in Fig. ^for com- 
paring fit results with the data for the calcium nucleus. The parametrization results 
and uncertainties are calculated at Q^=5 and 50 GeV^ for F2 and Drell-Yan, respec- 
tively. Because the data are taken at various points, the curves cannot be really 
compared with the data due to the difference. However, because the scaling vi- 
olation is not a large effect, it is obvious in Fig. ^that the analysis well reproduces 
the data. The F2 data play a role of determining antiquark distributions at small x 
and valence-quark distributions in the medium- and large- a; regions. On the other 
hand, the Drell-Yan data provide a constraint on the antiquark distributions at 
X - 0.1. 

From the analysis, we obtain the NPDFs for nuclei with A =2 to about 208. 
Among them, we show nuclear modifications of a medium size nucleus, calcium. 
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Figure 1. Fit results are compared with the data for F2°'/F^ and a^j^y / ^^dy ■ The curves and 
shaded areas indicate the fit results and their uncertainties at Q^=5 or 50 GeV^. 
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In Fig. [51 the weight functions are shown for the distributions u^", g*^", and g 
at Q^ = l GeV^. The shaded areas indicate uncertainties also at — 1 GeV^. The 
antiquark distributions at small x and the valence-quark distributions are well de- 
termined. However, the antiquark distributions at x > 0.2 and gluon distributions 
in the whole x region have large uncertainties. The valence-quark distributions are 
determined by the F2 data at medium x, and the small- a; part is constrained by 
the baryon- number and charge conservations. For fixing at large x and g^, we 
need new data which are sensitive to these distributions 5 . 
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Figure 2. Weight functions «>u„, wq, and Wg and their uncertainties are shown for the calcium 
nucleus at Q^=l GeV^. 



4 Effects on NuTeV sin^ 9w anomaly 

The NuTeV collaboration suggested that their neutrino scattering data should indi- 
cate a significant deviation of sin^ 9w value from other measurements. Their value 
is sin^ 9w — 0.2277 ± 0.0013 (stat)± 0.0009 (syst) in comparison with a global anal- 
ysis of other data: sin^ 9w = 0.2227 ± 0.0004. Because their experiments use the 
iron target, nuclear corrections have to be investigated. Among various nuclear cor- 
rections, it was shown in Ref. [3] that the difference between nuclear modifications 
of u„ and should contribute to the determination of siv?' 9w We investigate 
more details by using ^ analysis results for nuclear valence-quark distributions. 

For extracting sin^ 9w, the following Paschos-Wolfenstein (PW) relation is use- 
ful: R- = {a'^^^c - ^Nc)/i^CC - '^cc) = 1/2 - sin^ 9w, where and aX/^ indi- 
cate deep inelastic neutrino-nucleon cross sections for charged-current and neutral- 
current processes, respectively. This relation is derived for the isoscalar nucleon 
without strange-antistrange asymmetry (s — s — 0). We should note that the 
NuTeV experiments have been done for the iron nucleus which has a neutron ex- 
cess. There could be various nuclear corrections to the PW relation. In fact, writing 
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down the relation for a nuclear target, we obtain 



sm" t/iv ) {1 



ev{x) en{x)} 
1 4 



sin^ 0w {£vix) + £„(x)} 



sm Uw £. 



1 



fl 



— :^{ev{x) +en{x)] + 



3 

2{es{x) 



(1 



1 + £1,(0;) £n{x) 



l-(l-y) 

Here, the correction factors 



(4) 



£„, £s, and £c are defined by £„ = {wd^ — 

WuJ/iWd^ + WuJ, £« = ((A^ - Z)/A) {Uy - dy)/(Uy + Cty) , Eg = / [w y {Uy + Cly)] , 

and £c = c^/[wy [uy -\- dv)]^ where N is the neutron number and qy is defined by 
qv — q — Q- The functions Wu^ and Wd„ indicate nuclear modifications for Uy and dy 
distributions as explained in sectional The function Wy is defined by the average, 
Wy = {wd^ + 'Wu^)/2. We note that Eq. becomes the PW relation in the limit 

In particular, we investigate £„ effects on the b\v? 9w determination |3I4| . For 
calculating Sy effects, we should take the NuTeV kinematical into account. Such a 
kinematical effect can be incorporated by using functionals provided by the NuTeV 
collaboration. The nuclear modification (£„) effects on the sit? 9w are then calcu- 
lated. Our research is still in progress |4j; however, the preliminary result indicates 
a large uncertainty for the sin^ Ow determination. It means that the nuclear modi- 
fication could not be determined accurately and that there is still a possibility that 
the deviation is explained by the nuclear effect. 



5 Summary 

Nuclear parton distribution functions (NPDFs) have been obtained by analyzing the 
data for the structure function F2 and Drell-Yan cross sections. Their uncertainties 
were calculated by the Hessian method. The uncertainties are especially large 
for the antiquark distributions at a; > 0.2 and the gluon distributions, so that they 
need future experimental measurements. Such a determination of the NPDFs could 
contribute to the clarification of the NuTeV sin^ 9w anomaly. 
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